Geometric quantum gate for trapped ions based on optical dipole forces induced by 
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We present an implementation of quantum logic gates via internal state dependent displacements 
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I. INTRODUCTION 



Today, one of the most promising physical systems for 
an experimental realization of quantum computation is a 
string of trapped ions [jj] . In such a realization, a quan- 
tum bit (qubit) is represented by two internal states of an 
ion on which quantum logic operations can be performed 
through laser-ion interactions. Any quantum logic opera- 
tion can, e.g., be composed of single-qubit operations and 
ControUed-NOT (CNOT) gates between any two ions in 
the string J3] • While single-qubit gates are relatively sim- 
ple to perform, the CNOT and equivalent two-ion gates 
are more demanding and usually also slow as compared 
to single-ion gates. Only very recently such two-ion gates 
have been demonstrated experimentally SJi 9 ■ One class 
of two-ion gates is the geometric quantum gates, which 
can be relatively fast as compared to other two-ion gates 
and furthermore allows a high fidelity since the internal 
states of the ions are not directly involved in the gate 
operation Q,|3- 

In this paper we present a proposal for an experimental 
realization of a geometric quantum gate. The idea under- 
lying our gate proposal is that internal state dependent 
optical dipole forces applied to a pair of ions will displace 
their equilibrium positions in the trap in proportion to 
the forces acting on them. Their relative distance is thus 
modified, and the associated shift of their Coulomb po- 
tential energy gives rise to a complex phase factor which 
depends on the internal states of both ions, and which 
hence provides a non-trivial two-qubit gate. 

The paper is organized as follows. In Sec.^]we present 
a full quantum analysis of the dynamics of the system to 
obtain the correct values for the phase factors and to 
properly assess the effect on the atomic motion in the 
trap. In Sec. IIIII we consider specific implementations 
in alkaline earth ions. Assuming a Gaussian intensity 
profile of the applied laser beam, we find in Sec. llVl ex- 
perimentally relevant parameters for an implementation 
of the gate proposal using ''"Ca"'' or -'^'^^Ba^ ions and dis- 
cuss relevant error sources. Finally, in Sec. we give a 
discussion and a conclusion. 



II. THEORY 

We consider a pair of ions confined to the trap axis of 
a linear Paul trap, i.e., the axis where the rf- field van- 
ishes [3 ■ A laser beam which is far-off resonant with all 
internal transitions in the ions and propagates perpendic- 
ular to the trap axis induces an optical dipole force. This 
force can be described by a potential Udip{z, a) which de- 
pends on the position variable z along the trap axis, and 
on the internal state label a taking one of two possible 
values, represented by ] and [ in the following. The po- 
tential is controlled by the intensity, waist, polarization, 
and wavelength of the laser beam [3 ■ 

For an ion string illuminated by a dipole force inducing 
beam its total potential energy, U , is the sum of the 
dipole potential of each ion, the potential energy of the 
ions due to the trap and due to their mutual Coulomb 
repulsion. Specifically for two singly-charged ions of mass 
m, we have 



1 



U{zi,Z2,ai,a2) = -muj^{z^ + z^) + 



2 "'^■^''- ' '^' ■ 47reo|z2 - zx\ 

+ Vdi-p[z\^eq, ai) - Fdtp{zi^eq, ai){zi - Zl,eq) 
+ Udtp{z2,eq, 012) - Fdip{z2.eq, a2){z2 - Z2,eq), 



(1) 



where lUz is the single ion oscillation frequency in the 
trap and where a linear expansion of the dipole potential 
around the equilibrium positions Zi^eq (* — 1, 2) has been 
applied. Fdip{zi^eq,ai) = -dUdip{z,ai)/dz{z = Zi^^q) is 
the optical dipole force exerted on the i'th ion. 

Rather than considering a CNOT gate, we shall in 
the following focus on implementation of the equivalent 
ControUcd-Z two-qubit gate Q. A ControUed-Z gate 
applied to a superposition of the product states ||) ||), 
I i) IT), IT) I i) and IT) IT) changes the sign of the |T) |T) 
term (i.e., the phase of this state is changed by tt), while 
leaving the others unchanged. The implementation of a 
geometric Controlled-Z gate relies on the fact that the 
dipole force is internal state dependent. The origin of 
this dependence will be described in the next section. 
As mentioned above, the idea is that the dipole force 
displaces the ions away from their equilibrium positions, 
and the associated change in Coulomb energy gives rise 



to phase-shifts which depend on the internal states of the 
ions. By choosing a suitable temporal and spatial profile 
of the dipole force inducing beam, the obtained phase 
shifts of the above mentioned four product states can be 
made equivalent to a ControUcd-Z gate. 

To describe the motion of the ions, we first rewrite 
Eq. CO as follows 

U{zi^eq,Z2.,eq, Ckl, 02, t) = -mCJ^Az^ + -mcjj (z+ + 3zi) 

(2) 

+ Udip{zi^eq,ai,t) + Udip(z2,eq,Oi2,t) 
1 

y= [Fdip{z2,eq,Oi2,'t) — Pdip{zi^eq,Oii,t)] Z_, 



[Fdip{zi.eq, ai,t) + Fdip{z2,eq, "2, *)] Z^ 



where z+ = -75(^2 + 21) and z_ = -y={z2 — zi — Az) are 
the motional mode coordinates for the so-called center- 
of-mass mode and the breathing mode [Sj, respectively, 
and Az = Z2.eq — zi,eq dcuotcs the equilibrium distance 
between the ions. Turning to a quantum mechanical 
description, we introduce lowering (raising) operators a 
{a^} and b (b^) for the mode coordinates, i.e., z+ = 

y/h/{-imuJz){a + a'^) and z_ = ^Jh/{AmV3Lu^){b + b^), 
to find the Hamiltonian 



H = huj^a) a + vSfkUzb'b 

+ Udip{zi^eq,ai,t) + Udip{z2^eq,a2,t) 

+ f+{ai,a2,t){a + a^) + /_(ai, 02, i)(& + b'^), 
where 



(3) 



/+(ai,a2,i) = - 



8muj. 



and 



f-{ai,a2,t) = -. 



8m\/3u 



■[Fd2p{zi^eq,ai,t) (4) 

+ Fdip{z2.eq,a2,t)] 



-[Fdip{z2^eq,a.2,t) (5) 

— Fdip{zi^eq,ai-,t)] 



are responsible for excitation of the center-of-mass mode 
and the breathing mode, respectively. The time- 
evolution of the system can be described by a uni- 
tary time-evolution operator hi, which evolves state vec- 
tors * in time according to ^(t) = U{t)^{t = 0). 
U{t) solves the time dependent Schrodinger equation, 
im{t) = H{t)U{t), U{0) = /. Since the Hamiltonian 
can be decomposed as a sum of commuting terms, lA can 
be expressed as a product 



U =exp 



xU+(t)U-{t) 



dt' [Udip{zi^eq,ai,t') + Udip{z2,eq,a2,t')] 

(6) 



where U+ (t) and U- (i) are time-evolution operators cor- 
responding to the Hamiltonians 



and 



H^ =hLUzCi^a + f+{ai,0i.2,t){a + 0,^) (7) 



H_=V3hw,b^b + f^{ai,a2,t){b + b^), (8) 



respectively. In the following we solve for U+ with the 
understanding that the solution U- can be obtained from 
14+ simply by replacing /+ with /_ and uj^ with \/ouJz- To 
this end, we switch to the interaction picture with respect 
to the Hamiltonian of the free harmonic oscillator 






(9) 



and make the Ansatz that 

ZY™t,+ =6*"^*"^"^^+ = e-l'5+l'/2e#+e^/3;«*e*'3+''. (10) 

From the Schrodinger equation for Uint,+ it follows 
that the harmonic oscillator phase-space displacement 
(3+{ai,a2,t) = p+fy/mfkOz - iz+/ yjh/{mujz), with p+ 
being the center-of-mass mode momentum, can be writ- 
ten as 

1 /■* 
l3+{ai,a2,t)=-- dt7+(ai,a2,t')e"'"'' (H) 

and that the phase (p+ acquired due to excitation by the 
force term /+ is given by 



6+(ai,a2,i) =- -p- Im 



dt' f+{ai,a2,t')e 



f\ ^ — iu)zt' 



(12) 



X / dt"f+{ai,a2,t")e 



ir\iLj-t" 



The displacement of the ions is internal state depen- 
dent, which leads to coupling (or entanglement) between 
the internal and the motional states. For the gate op- 
eration, this is an undesired effect and we shall there- 
fore request the displacement to be zero at the end of 
the gate operation, such that U, 



int,+ 



e"P+. In the 

implementation of the ControUed-Z gate described be- 
low, we will take the dipole-potential to be of the form 
Udrp{zi,ai,t) = Uconstizi) + Uosc{zi,ai)g{t), {i = 1,2) 
in a time interval [0, T] and zero otherwise, specifically 
with g{t) = outside this interval. To avoid the large 
internal state dependent phase factors in Eq. © we fur- 
thermore assume /„ g{t)dt = 0. The force term /+ 
can be written as /+(ai,a2,i) = fconst{zi,eq-,Z2,eq) + 
fosc{zi,eq-, Z2,eq, ai,a2)g(t), which together with Eq. ^^ 



implies that the center-of-mass mode displacement at the 
end of the gate operation can be written 



fosc [^ 

i'Y' / dtg{t)e' 



nuJr. 

(13) 

Both terms of this expression vanish if T is an integer 
number n of oscillation periods, T = I'nnjiOz-, and if the 
Fourier transform 



){^z) = 



1 



/27r 



dtg{t)e 



— iuJzt 



1 



dtg(t)e 



(14) 



= 0. 



From Eq. H12|) it can be shown that the phase acquired 
during the gate operation can be expressed as 



MT)=^Jd.'-^ 



/fK) 



= C'l/co«.t + C2fconstIoscg{^) + T2 ^UJ — 



fi2 



(15) 



UJ — UJz 



where /+(w) is the Fourier transform of f+{t) and Ci 
and C2 are constants. The first term on the r.h.s. is 
irrelevant, since it does not depend on the internal state 
of the ions and the second term vanishes since we have 
required 5(0) = {2Tr)~^^'^ J^ g{t)dt = 0. The interesting 
term is the last one, from which we observe that the closer 
the characteristic frequencies of the function g{t) are to 
the oscillation frequency coz (or vS'^z for the breathing 
mode), the larger is the accumulated phase. It is thus 
a natural choice to adopt a harmonic time dependent 
force that oscillates with just one oscillation cycle less (or 
more) than the trapping motion during the interaction 
time T. In Fig.QIa), one example of a phase-space trajec- 
tory of the center-of-mass mode of two ions in the ||) ||) 
state is shown. The duration of the interaction is taken 
to be 15 trapping periods, and 14 periods of the applied 
periodic force. We observe that the net displacement 
vanishes at the end of the gate operation. For the same 
interaction, Fig.^b) shows the phase-space trajectory of 
the breathing mode with the ions being in the |i)|t) state. 
The breathing mode frequency of V^uJz is far-off reso- 
nant with the frequency of the applied force and hence 
the breathing mode is much less excited than the center- 
of-mass mode. Note that since 15\/3 — 25.98 ^ 26, this 
oscillator mode undergoes an almost integer number of 
oscillations and hence it experiences a nearly vanishing 
net displacement at the end of the gate operation. Other 
good choices for the duration of the interaction in units 
of the trapping period are 56 (56-\/3 = 96.995 ~ 97), 
and 209 (209^/3 = 361.9986 ~ 362). The phases 0± are 
equal to the areas of the {z±,p±) phase-space trajecto- 
ries in units of h, and they depend on the internal state 
of both ions, as needed for a two-qubit gate. Fig. ^c) 
shows the build-up of the phase during the interaction. 



The parameters have been chosen to ensure an effective 
phase shift of n on the |t) |T) state (see Eq. (I55|) below) as 
desired for the Controlled-Z gate. In the following section 
on the physical implementation of the scheme, analytical 
expressions for the acquired phases are provided. 



III. IMPLEMENTATION IN ALKALINE EARTH 

IONS 

A. Atom-light coupling 

Since alkaline earth ions are among the most promi- 
nent atomic ions for doing quantum logic, we will now 
focus on experimental realizations of the Controlled-Z 
gate based on such ions with presentations of specific re- 
sults for ^°Ca^ and ^"^^Ba"''. We choose the qubit states 
to be li) = n^S,/2i-l/2) and |t) = n^S./^i+l/^), 
on which single qubit operations can be performed us- 
ing Raman transitions via the n'^Pi/2 level. Assuming 
that the frequency lol of the dipole force inducing laser 
is close to or below the transition frequencies LO112 and 

(^3/2 of the n'^Si/2-n'^Pi/2 and the n'^Si/2^n'^P-i/2 tran- 
sitions, respectively, we may to a good approximation 
only consider contributions from these two transitions to 
the dipole potential (as long as ojl is not in the imme- 
diate vicinity of the transition frequencies of the weak 
'T'^'S'i/2 ~ (^i ~ 1)^^3/2,5/2 electric quadrupole transi- 
tions). Expanding the dipole force inducing beam into 
cr+- and (T~-polarized light components with respect to 
its propagation axis (the quantization axis), the contri- 
butions to the dipole potential for the two states ||) and 
\]) arc those associated with the transitions indicated in 
Fig. 121 The respective dipole potentials can consequently 
be written 

Ui = ^+1+ + ip-I- and [/| = i'^I+ + tP+I^ , (16) 

where I± is the intensity of the a'^- and cr^-polarized 
components, respectively, and where 



^^ 



3nc^ 



2r 



1/2 



3uj^ 



1 



1 



/2 \^l/2-~ ^L ^1/2 +^L 



(17) 



- 3/2 



3uji 



3/2 \^3/2-^L ^3/2+W-L/ 



and 



^- = 



37rc2r 



3/2 



=^3/2 



1 



1 



^3/2 ~ '^L W3/2 + ^L 



(18) 



depend only on the properties of the ion and the fre- 
quency ujl of the dipole force inducing laser. Here, ri/2 
and r3/2 are the transition strengths of the n'^Si/2- 
n'^Pi/2 and the n'^Si/2-n'^P3,/2 transitions, respectively. 
As can immediately be seen, the force derived from this 
dipole potential will be different for the two qubit states. 
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FIG. 2: Relevant energy-levels and transitions in alkaline 
earth ions (e.g.,*°Ca+, ®*Sr+, and ^^®Ba+) for calculating the 
dipole potential for the qubit states |J,) — ^Si/2(— 1/2) (solid 
lines) and \X) — ■^ 51/2 (-1-1/2) (dashed lines) due to a dipole 
force inducing laser beam containing cr"'"- and cr~-polarized 
components with intensities /+ and /_ , respectively. ■(/)+ and 
'0- are defined in the text. 



whenever the intensities of the two polarization compo- 
nents differ. Hence, by introducing a sinusoidal temporal 
variation of the intensity of the polarization components 
given by 

/±(z,i) = ^/(;^)[l±sin(m)], (19) 

which can be done by, e.g., using an electro optic phase 
modulator, the situation considered in Sec. ^ can be ob- 
tained. 

From Sec. m we thus identify g[t) — sin(rit) for < 
t <T. As for the choice oi fl, we recall the requirements 
T = 2TTn/u!z with n integer, g{0) = g{^z) — giV^^z) — 
and that the characteristic frequency of g{t), i.e., 17, 
should be close to ujz- Consequently, a good choice is 
i7 = (1 — l/n)u!z with n 3> 1, as in the example of Fig.Q] 
where n — 15. The Fourier transform of g(t) on the 
interval [0, T] contains two terms which are proportional 
to sin[(w - n)T/2]/{u; - n) and sin[(w + n)T/2]/{u; + n), 
respectively. These sine-functions peak at w = O and 
Lu — —fl, while they have exact zeros at lu — and lu ~ 
uJz- Furthermore, they are suppressed a,t uo — VSij-'z, 
especially when VS't- is close to an integer. 



FIG. 1: Real and imaginary parts of the center-of-mass mode 
and the breathing mode displacement /3± (f ) and the acquired 
phase with n = 15. The parameters are chosen such that 
an effective phase shift equal to tt of the |t) |t) state (see 
Eq. 1331 below) is obtained. (a) Parametric plot show- 
ing (Re[/3+(f)],Im[/3+(t)]) = {p+/^/hmujz, -zj,I ^hlmi^:,) 
when the ions are in |i)|i). (b) Parametric plot showing 

(Re[/3_(t)],Im[/3_(i)]) = {p- j \J hm^fZ<jj,.~z- j sjhlm^/?,u.?, 
when the ions are in |i) |T)- In both (a) and (b) the phase- 
space trajectory starts out from the origin and goes clockwise 
as time elapses, (c) Acquired phase. The solid line is the to- 
tal effective phase of 2[^+(j,J,, t) — </)_(||, t)\. The dashed line 
is the center-of-mass mode contribution of 2(/)+(J,J,, i). The 
dotted lines indicate phases of and tt. 



B. Phase shifts 

In order to calculate the displacement and the ac- 
quired phase for the center-of-mass mode and the breath- 
ing mode, respectively, for the specific intensity variation 
in Eq. I|I9|) . the /-functions entering in the integrals in 
Eqs. (|ll|l and (|12|) have to be determined using the defi- 
nitions in Eq. Q or ((S)). For /_(-, which is responsible for 
excitation of the center-of-mass mode, we find 



/+ai)=-[/o+ + /i+sin(m)], 
/+(iT)=-[,/o+ + /2+sin(m)], 
/+(Ti)=-[/o+-/2+sin(m)], 
/+(tT)=-[/o-f-/i+sin(Oi)], 



(20) 
(21) 
(22) 
(23) 



where 



f2. 



8muj: 



■(Fi+F2)(V++^-), 



Smuj- 



■(Fi+F2)(V+-^-), 



-(Fl-F2)(V^+-^-) 



and 



F,.^- 



Idliz) 
2 dz 



(z = l,2). 



(24) 
(25) 
(26) 

(27) 



The full time-dependent expressions /3(i) and (pit), which 
were used for obtaining the plots in Fig. ^ can be found 
by carrying out the integrals in Eqs. (|ll|l and p2|l for the 
center-of-mass mode as well as the breathing mode [lOJ. 
Here we only state the acquired phases at the end of the 
gate operation, which for the center-of-mass mode are 
given by 






fl 



2uj,T 



Pu 



(2/1)2 ^2 _ q:2 (;j^^)2 



(28) 



/2+ 



2w^r 



/!+ n^^ 



(2fi)2 a;2 - f]2 Q^^y 



(29) 



neglecting a term which is independent of the internal 
state and assuming n ^ 1. These phases scale quadrat- 
ically with n with one factor of n originating from the 
time T and the other originating from the denominator, 
due to the fact that Q. was chosen to be near-resonant 
with LOz- The phases acquired due to excitation of the 
breathing mode are likewise found to be 



</)„(ii,r)=^_(tT,T) = 
<^_(iT,r)=0_(ti,T)- 



/f_ 2^/iwzT ft VSTi 



mr 


3^2- 


r!2 ■ 


" {fy^z? 
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(30) 


fL 


2V3t. 


zT 


. /I- 


VonTT 



where 



h- - /i 



F2 — Fi 

v^{F2+Fi) 



(2/1)2 3^2 _ Q2 



and /2- = /: 



(r^^)2 2 

(31) 



i^2 + -Fi 



2+ 



^(^2 



^1) 

(32) 

These phases only scale linearly with n, since the rota- 
tion frequency Q is off-resonant with the breathing mode 
frequency of \/^uJz- In Eqs. H30|l and H31|l. terms which 
are smaller than the stated terms by a factor of n and 
further suppressed if VSn is close to an integer have been 
neglected. 

The combined effect of the above phase shifts is equiv- 
alent to a single phase shift (f>± = 4'±{ii) — 0±(iT) — 



Mn) + Mm = 2[0±ai) - ^±(1?)] of the mir) 

state [III. Thus, to make a ControUed-Z gate, we re- 



11+. 

m 

quire that 



^=2[</)+ai,T) + 0_(ii,r)-<^H.(it,T)-(/._(it,T)] 

(33) 

/l+ ^ /2+ + ~~~ (A- ^ /2- 



{fllOz 



IV. REALIZATION OF THE CONTROLLED-Z 
GATE USING GAUSSIAN LASER BEAMS 

A. Intensity requirements and ofF-resonant 
scattering rates 

In the following we will consider possible realizations 
of a ControUed-Z gate in a two ion string using one or 
two dipolc force inducing laser beams propagating per- 
pendicular to the trap axis as sketched in Fig. |31 In the 
first case shown in Fig.|2Ia), we assume that the equilib- 
rium distance Az between the two ions is smaller than 
the waist of the laser beam, such that they feel essentially 
the same strong dipole force. The situation in Fig. ^h) 
corresponds to a case where the size of the force on the 
two ions are the same but of opposite sign for identi- 
cal internal states. Finally, Fig. ^c) shows a realization 
similar to the first case, but with application of two laser 
beams. 

Since in all cases we will assume the ions to be situated 
close to one of the points in the intensity profile where 
the induced dipole force is largest, much of the analysis is 
similar for the three situations. Hence, in the following 
we will focus on the situation in Fig. Ola) in order to 
establish formal equations. Assuming the ion string is 
centered at z = 0, we may write the intensity of the 
dipole force inducing beam as 



I{z) = loe 



-2iz-zoy/W^ 



(34) 



where zq = W/2 is the center of the dipole force inducing 
beam, W is the beam waist and /q is the peak intensity. 
For the gate operation two highly relevant and linked 
parameters are the required peak intensity Iq and the 
fidelity loss due to off-resonant scattering events. The 
required peak intensity enters in Eq. l|55)l through the /- 
functions. For a given wavelength and a given choice 
of beam waist the required peak intensity can be de- 
termined. Using the definitions in Eqs. l(^ . l(^ . l(?7|) 
and (Eg, we find i?i « F2 « e'^/'^Io/W, /2+ « /i_ « 
and a pair of non-zero expressions for /i_|_ and /2- , which 
together with Eq. H33|) yields the required peak intensity 



/o 



/ 2e^fiujlmW^ 



(35) 



where only the leading term in n has been retained. The 
scaling with the various parameters is intuitively reason- 
able. First, it has already been discussed that large n 
leads to a large phase pick-up and hence a low intensity 
requirement. Second, the larger the difference ip^ — tp^, 
the larger is the dipole force difference for |j) and \1). 
Third, if the waist is small, the dipole force is large, which 
in turn reduces the required intensity. Note also that the 
required laser power (~ IqW^) is proportional to W^, 
which makes a small waist very attractive. 

Knowing the required intensity, we can now determine 
the average probability Pgc for a scattering event from 
one of the two ions during the gate operation. Formally, 
we can write Psc ~ ^scT {T^cT <C 1), where Tsc is the 
average total scattering rate of the two ions. Assuming an 
equal average population of the two internal states and 
an intensity of e^^'^Io at the position of the ions, the 
dominant scattering from the n^Pi/2- and the n^P^/2- 
state yields 



(a) 



P.. 



where 
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(36) 



(37) 



The dependency of P^c on the internal structure of the 
ion and the wavelength of the dipole force inducing laser 
is contained in the front factor of Tsc/{''P+ — ^-), show- 
ing that Psc can be minimized either by making Tgc 
small or ^+ — tp^ large. T^c becomes small in the limit 
ojl ^ Wi/2,^3/2 due to the factor of a;£, however, in 
the same limit ■0+ — V'- is proportional to the difference 
ri/2^1/2 ~ ^3/2^3/2' which is also small. Alternatively, 
if uji/2 < lul < ^^^3/2: 4^+ ~ "0- can become a sum of 
two positive terms, however, for Fjc to be small in this 
case, a large fine-structure splitting is required. Both in 
the far-off resonant case and when the dipole force in- 
ducing laser is tuned in between the fine-structure levels, 
the ^'^^Ba^ ion turns out to be more attractive than the 
'^''Ca^ ion. Another candidate ion is ^^Sr"*" which for the 
present gate proposal is less attractive than ^■^'^Ba^ but 
more attractive than ^''Ca , however, in the following we 
only consider ''"Ca'*' and ^^^Ba^. 



B. Experimental parameters for ^''Ca^ and ^"^^Ba^ 

To get more quantitative numbers out, we now consider 
the cases of two '^''Ca^ or two ^'^^Ba'*' ions in a trap with 
ti;^ = 27r X 1 MHz, which leads to an equilibrium distance 
Az of 5.6 /im and 3.7 /xm between two ions, respectively. 
First we consider the situation depicted in Fig. Ola) for 
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FIG. 3: Three possible configurations for realizing the two- 
ion gate considered in the text, (a) Two ions positioned 
at distances of W 12 ± /S.z/2 from the center of a Gaus- 
sian laser beam, (b) Two ions positioned at distances of 
±W/2 — ±Az/2 from the center of a Gaussian laser beam. 
(c) Two ions, each positioned at a distance of W/2 from the 
center of a tightly focussed Gaussian laser beam. 



Az ^ W, such that the dipole forces on the two ions are 
equal. To fulfill this we choose a waist of 30 /xm. Finally, 
by choosing n = 15 (von « 26) the intensity and the 
scattering rate can be calculated for a given wavelength 
of the dipole force inducing laser beam. For ^°Ca~'' and 
^^^Ba^ with the dipole force inducing laser tuned either 
in between the fine-structure levels or far red detuned, 
we find the following values for (Power, Pg^, A^), where 
Xl is the dipole force inducing laser wavelength. '^'^Ca^: 
(8W, 30%, 395.1 nm) and (~0.5MW, < 4%, >1500nm). 
"8Ba+: (86W, 6%, 474.5nm) and (-33kW, < 1.2%, 
>1000nm). Clearly, this is not very promising for exper- 
imental realizations. However, since the main problem 
is the requirement that W 3> Az, the situation depicted 
in Fig. Ofb) is more favorable. Here W = Az with the 
ions positioned symmetrically around the center of the 
dipole force inducing beam. Since we do not have the 
requirement that Az ^ W^ we can choose a lower trap 
frequency of coz — 2tt x 200 kHz, which implies that Az 
is equal to 16.4 /urn for "^^Ca^ and 10.9 /im for -'^•^^Ba'*' 
and which yields the following values for (Power, Pgc, 
Xl). ''°Ca+: (120 mW, 8%, 395.1 nm) and (-6.5 kW, 
< 0.9%, >1500nm). 1383^+: (360 mW, 1%, 474.5 nm) 
and (-140 W, < 0.2%, >1000nm). For high fidelity gates 
the required power is experimentally still too demanding. 
Finally, we consider the situation presented in Fig.|2Ic) 
where the waist of the applied laser beams is not directly 
related to the equilibrium distance between the ions and 
hence rather tightly focussed beams can be applied. As 
long as the beam waists are larger than the thermal ex- 
cursion of the ions from their equilibrium positions, the 
equations derived above for the case of Fig.|3Ja) is imme- 
diately applicable. For uiz — 27rx 200 kHz and W — 5 fim, 
we plot in Fig.^the required laser power and TscT ~ Psc 



for *°Ca~'' and ^'^^Ba''' as a function of the dipole force 
inducing laser wavelength. The plots in Fig. 01a,c) ex- 
tend to a wavelength of 5 fim in order to show the long 
wavelength behavior, however, one should keep in mind 
that due to diffraction it is technically demanding to ob- 
tain a 5 /im waist for the longest wavelengths. In all plots 
n = 15 was chosen, which gives a very reasonable gate 
time of T = 75 /zs. 



For ^°Ca''', we see from Fig. 0fa) that a high power of 
~ 200 W is required in the long wavelength limit for the 
considered parameters. However, by choosing n = 209 
instead of 15, the required power drops to '--^ 15 W at the 
cost of an increased gate time of T = 1 ms. Considering 
the challenge of diffraction a laser at a wavelength of 2 /im 
would be favorable in the long wavelength limit. Fortu- 
nately, the required power of ~ 15 W is easily achievable 
with commercial Thulium fiber lasers operating in the 
range of 1.75 - 2.2 /xm [l^. Hence, high fidelity gates with 
Psc ~ 10"'^ should be possible. Other interesting lasers in 
the far-off resonant regime are the Nd:YAG (yttrium alu- 
minium garnet) laser at 1064 nm and the CO-laser near 
5 /im, but the relatively high scattering probability and 
diffraction limitations, respectively, makes these lasers 
less attractive than the Thulium fiber laser. A wave- 
length near 395 nm is furthermore attractive for ^°Ca''' 
(see Fig. 0Jb)), since the needed power of only ^ 3mW 
easily can be obtained. The scattering probability is a 
few percent, which is acceptable for a first demonstra- 
tion, but not good enough for implementation of error 
correcting schemes. By a further reduction of the beam 
waist, e.g., by placing a lens system inside the vacuum 
chamber where the ion trap is situated, it should be pos- 
sible to focus to below 1 /im [l3| and hence reduce the 
scattering rate by a factor of 5 — 10. 



Due to the larger fine-structure splitting of ^"^^Ba^, 
the scattering probability is already below 1% for a 
wavelength around 475 nm (see Fig. 0Jd)) between the 
Q^Si/2-Q^Pi/2, ^-p3/2 transition wavengths. The low 
laser power needed in this wavelength region can easily be 
supplied by frequency doubled diode laser systems. Note 
also that light from an Argon-ion laser at a wavelength of 
488 nm could be a reasonable possibility. However, as for 
^^Ca"*" a Thulium fiber laser seems to be most ideal, since 
high fidelity gates with Pgc ~ 10~^ (a value comparable 
to the threshold value required for fault-tolerant quan- 
tum computation Il3l | ) should be feasible with laser pow- 
ers of about 14 W for the case with n = 15 and T ~ 75 fj,s 
as shown in Fig.QJc). In the long wavelength limit, the 
1064 nm wavelength of a Nd:YAG laser could also be a 
possibility since a scattering probability of ^ 10^"^ at a 
laser power of ~ 7W is achievable. It should be noted 
that in all cases discussed above, the required power, the 
scattering probability and the gate time can be adjusted 
by changing n, ujz and W. 



C. Error sources 

When we considered the realization of the gate schemes 
above, we already discussed the effect of spontaneous 
emission on the fidelity of the gate operations. There are, 
however, other experimental error sources which will also 
influence the final fidelity of an actual realization. A key 
element in the present gate proposal is the polarization 
rotation, which serves to remove phases due to internal 
state dependent Stark shifts. Since the Stark shifts are 
of first order in intensity and the phase acquired by the 
change in Coulomb energy only originates from second 
order effects, even small errors in the laser parameters 
may be very critical to the actual gate operation. In 
the following subsections, errors due to non-perfect bal- 
ancing of the two polarization components of the dipole 
force inducing laser beam, timing errors as well as power, 
position and frequency- fluctuations of the dipole force in- 
ducing laser are discussed. 



1. Polarization errors 

In case there is an imbalance of the intensity in the 
two polarization components, such that 

I± {z,t) = h[z) [1 ± ep] [1 ± sin[nt)] , (38) 

where ep accounts for the imbalance, there will be two ex- 
tra terms in the dipole potential. One term [oc ep sin(r2t)] 
enters with the same sign in C/| and U^ and hence it does 
not give rise to any Stark-shift induced phase-difference 
between ||) and jt). The other term, epI(z){ip+ — tp-)/2, 
enters with a different sign in U^ and C/| , which leads to 
a phase-difference of Acj) = epl{z){^+ — ■i/j^)T/h between 
II) and It) in each ion at the end of the gate operation. 
Since this difference should be much smaller than the de- 
sired phase-shift of tt, we find by using the expression for 
A(f), Eqs. 1(231), ^B and T = 2-1:11 /uj^ the condition 



e«< 



SujzmW^ 



(39) 



In the situation considered above with W = 5 /tm and 
LUz ^ 2t: X 200 kHz this means ep < 1/400 for '*°Ca+ 
and ep < 1/700 for 1283^+. FulfiUing these criteria 
is not very easy, but fortunately, the undesired phase- 
difference can be cancelled by using a type of spin-echo 
technique /14|. Instead of generating the full effective 
phase-shift of tt in a single operation, the gate-operation 
can be performed in four steps: (1) Run the gate at half 
the intensity, to get an effective phase-shift of 7r/2 on 
It) It)- (2) Swap the population between ||) and |t) by 
applying single-qubit 7r-pulses to both ions. (3) Same as 
(1). (4) Same as (2). The trick here is that the unde- 
sired phase-differences due to a polarization error, which 
are obtained in step (1) and (3), are of the same mag- 
nitude but have opposite signs (due to the population 
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FIG. 4: Required power and TgcT vs. wavelength of the dipole force inducing laser for '*''Ca+ and '^^^Ba^. In all plots 
LUz — 2tv X 200kHz, W — 5/iin, n = 15 and T — 75 pis. (a) ^"Ca^, far red detuned laser. The divergences near 900 nm 
are due to a cancellation of tp+ — ?/>_. (b) ^^Ca"'', laser wavelength in the vicinity of the 4 '^ 51/2-4 ■^Pi/2 and 4 ^51/2-4 ^P3/2 
transition wavelengths, which are indicated by the symbols ^Pi/2 and ■^P3/2. (c) ^'^^Ba"'", far red detuned laser, (d) ^^®Ba+, 
laser wavelength in the vicinity of the 6^5*1/2-6^^1/2 and 6'^5i/2-6'^P3/2 transition wavelengths. 



swapping) and therefore cancel out. The gate operations 
in (1) and (3) both give an effective phase-shift of tt/2, 
even though the population is swapped in (2), because 
^±ilLT) = 0±(TT,r) and 0±aT,T) = 0±(n,T). The 
final TT-pulse just swaps the population back. Using this 
trick the gate time is doubled (neglecting the duration of 
the relatively fast 7r-pulses), while the required intensity 
is halved. Since the required intensity is proportional to 
n~^ and the gate time T oa n, these parameters can be 
re- adjusted if an appropriate value for n is available. 

Finally, the imbalance also gives rise to errors in the 
gate operation, which the spin-echo trick does not cancel. 
These errors are of order ei or ep/n and hence they are 
suppressed to the 10"'' level at a small but realistic value 



of. 



1%. 



2. Timing errors 

In case the gate time differs from the duration of a 
full number of polarization rotation periods, an unde- 
sired phase-difference will again build up. Assuming 
T = 6T + 27r(n - l)/rj and ^T < 1, the phase- 
difference is equal to A^ above, with the replacement 
ep i-> 5T'^Vt/{2T). Again using W — 5 fj,in and Uz = 
27r X 200 kHz and considering '^"Ca+ with n = 15. Eg. ifS^ 



translates to ST <C 0.5 /xs, with a more relaxed limit for 
larger n (the limit is proportional to ^/n). Applying elec- 
tro optic modulators to control the laser pulse length, the 
condition on 5T is not very severe since switching times 
of a few nanoseconds can be obtained. 



3. Power fluctuations 

If the total laser power fluctuates at a frequency LUf, 
such that the total intensity is given by Io{t) = /o[l + 
ef sin^LUft)], then the resulting fluctuations in the dipole 
potential integrated over the gate time T will give rise to 
a phase difference between ||) and |t). If the intensity 
fluctuations are random, they can in general not be ex- 
pected to cancel using the spin-echo trick. When tof ^ fl 
the phase difference Acj) ~ efl{z){4>+ — ip^)/2, i.e., the 
same phase-difference as above, just with ep replaced by 
£/, which means that e/ <C 1/400 is required for ^"Ca"''. 
When ujf ^ fl the phase-difference is smaller by a factor 
of 27m and even smaller if tuj ^ f2. Intensity stabiliza- 
tion fulflUing €f <C 1/400 is not unrealistic, in fact a com- 
mercially available Laser Power Controller already offers 
a power-stability of 3 • 10~^ within certain limits [13. 



4- Position fluctuations 

Fluctuations in the position of the dipole force induc- 
ing beam give rise to intensity fluctuations, which lead 
to imperfect cancellation of the Stark-shift induced phase 
and an unwanted variation in the dipole force exerted on 
the ions. Hence, for high quality gate operations with 
errors at the 10^** level, position jitter has to be of the 
order of 1% of the beam waist W or smaller. This condi- 
tion is most restrictive for the situation shown in Fig|3Jc), 
where a pointing stability of about 50 nm is needed when 
W = 5 /im. Though small, this type of stability is tech- 
nically possible. 



5. Frequency fluctuations 

As for the power and position fluctuations, laser fre- 
quency fluctuations will lead to fluctuations in the dipole 
potential and hence to an imperfect cancellation of the 
Stark-shift induced phase. Since laser frequencies can 
be very accurately controlled and frequency fluctuations 
anyway are expected to be small, as compared to the de- 
tuning from any of the two fine-structure levels, this is 
not expected to play any significant role. 



V. DISCUSSION 

The gate proposal presented above has some similar- 
ities with the gate recently demonstrated by the NIST 
group 4]; in fact the mechanism which gives rise to 
the desired phase shift is exactly the same. There are, 
however, also some essential differences between the two 
schemes. 

In the NIST experiment, a dipole force along a given 
trap axis was provided by the intensity-gradient of a 
moving standing wave light-field. The wavelength of the 
light-field and hence the period of the standing wave is 
set by the requirement that the Stark-shift induced phase 
shift of the two qubit levels should be zero (corresponding 
to ip+ + ip- = in our case), which is fulfilled if the laser 
is tuned in between two fine-structure levels. This can, 
however, give rise to a significant amount of scattering 
events, as we also saw above. Even though this scheme 
relies on a moving standing wave, the instantaneous force 
on the two ions is required to be equal. Consequently, the 
ions have to be well localized which requires cooling to 
the Lamb-Dicke limit with respect to the wavelength of 



the light-field, which in this case also is the Lamb-Dicke 
limit for the qubit operations. An equal dipole force on 
the two ions was obtained by adjusting the distance be- 
tween them to an integer number of standing wave pe- 
riods, which may be difficult to generalize to perform a 
gate between any two ions in a multi-ion string. 

In the proposal presented here, where the dipole force 
is provided by a variation in the beam-profile, the excur- 
sion of the ions from their equilibrium position should 
only be smaller than the beam waist, which is adjustable, 
but typically up to ten times larger than a relevant transi- 
tion wavelength. This means that except for very tightly 
focussed beams the Lamb-Dicke limit criterion need not 
be fulfilled and, specifically, the ions need not be cooled 
to the motional ground state. Since the dipole force in- 
ducing beam propagates perpendicular to the ion-string 
in our proposal, addressing of specific ions for implemen- 
tation of gates in a multi-ion string should be feasible. 
A theoretical description of this situation should also 
be quite straightforward by generalizing Eq. ©. Fi- 
nally, owing to the polarization rotation method, the 
dipole force inducing laser is allowed to be far-off res- 
onant with respect to the relevant internal transitions, 
such that a scattering probability below the asymptotic 
threshold value required for fault-tolerant quantum com- 
putation p^ in principle can be obtained in the long 
wavelength limit. 

It should be mentioned, that the idea of using optical 
dipole forces for implementing a ControUed-Z gate has 
also been considered by Sasura and Steane for an a rray 
of microscopic ion traps with a single ion in each trap [ll| . 
In many respects this system is very similar to the one 
considered here and naturally many of the considerations 
are the same as those made above. 

In a very recent proposal, Garcia-RipoU, ZoUer and 
Cirac present a geometric gate, where the momenta of 
the involved ions are controlled by absorption of photons 
from a discrete set of laser pulses [l^. In this case re- 
quirements on the pulses naturally arise for having zero 
displacement and for obtaining the desired phase shift. 
It can be shown that these requirements are discrete ver- 
sions of those expressed through g{io) above. 

In conclusion, we have presented a proposal for a ge- 
ometric quantum gate and shown that an experimental 
realization using the attractive alkaline earth ions *°Ca^ 
and ^^^Ba''" is feasible with gate times below 100 /iS and 
errors at the lO^'* level as required for fault-tolerant 
quantum computation |l3l | . 
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